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ANALYTICAL AND NUMERICAL. INVESTIGATION OF STRUCTURAL 
RESPONSE OF COMPLIANT 'WALL MATERIALS 


By 

R. Balasubramanian^ 


SUMMARY 


Surface motion of compliant walls in. drag reduction experiments has been 
analyzed. Critical comparison is. made between the dynamic motion of the structure 
and the postulated mechanism of drag reduction (ref. 1) . The spectrum of surface 
motion indicates that membranes over deep' cavities respond at low frequencies 
and large wavelengths. The membrane over a deep cavity is therefore found not 
to yield the desired response predicted by the postulated -mechanism. The mem- 
brane over a thin air gap is found to. act -as a wavelength chopper, and analysis 
of the nonlinear response of that compliant surface indicates its possible 
suitability for compliant wall experiments ► Periodic structures are .-found to 
lock in the desired wavelengths of motion, and it is found that at. least in 
Kramer’s (ref. 2) initial experiments they could have produced high frequency 
surface motions. Laminated structures- are found to be vexy ineffective as com- 
pliant models, except when there is no bonding between the membrane and the 
backing. Computer programs developed for these analyses, are documented in this 
report. 


INTRODUCTION 

The impetus for the compliant wall drag reduction program at -NASAr Langley 
Research Center, came mainly from, the. fact.-that -reduction in-turbul:ent..skin friction 
drag reported in the literature (refs. 2 to 5) may be translated into potentially 
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large savings in energy for CTOL aircraft.’ The group effort- at Langley was 
initiated with the goals of: 

1. Understanding the conditions ^ if .any, for favorable compliant motions, 

2. Diq)licating previous data under, rigorously controlled test conditions, 
and 

3. Designing improved experiments to further understanding of the compliant 
wall drag reduction phenomenon. 

The present work under NSG 1236 eminated from a desire to evaluate the struc 
tural response of compliant walls in successful drag reduction experiments. 

THE' PROPOSED FLUID -STRUCTURE INTERACTION MODEL 

Bushnell (ref. 1) has. proposed a compliant .wall interaction model by xvhich 
a fully turbulent flow over a vibrating compliant wall could produce reduced skin 
friction on the wall. In the past, few years, researchers have developed an under 
standing of some details of the streaky flow adjacent to. a rigid wall. The pro- 
posed mechanism of Bushnell (ref, Ij suggests that the induced aerodynamic 
pressure field due to the structural motion could impose a high frequency stab- 
ilizing modulation of the "preburst", flow and thereby reducing the rate of pro- 
duction of turbulent shear. Numerical e^eriments by Orszag (ref. 6) and active 
wall experiments reported by Kendall and. Collins (ref. 7j suggest that such a 
mechanism may well be valid.’ 

Based on this model, the favorable structural motion has a wavelength on 
the order of X* = 100, in wall units, [X = X* v/ (,i/C^/2 U_^j; where Cp is’ the 
skin friction coefficient, v is the kinematic- viscosity, and is the mean 

flow velocity], a wave speed C = 0..3 U^, and an amplitude a ss 0(1,0). For low 
speed air experiments, this criterion indicates that the drag reduction effect 
can only be observed in compliant walls -made of elastomers and carefully designed 
to obtain the desired wavemotions under given flow conditions.. Indeed, most 
of the reported "successful" experiments were in. compliant walls using thin 
membranes . 

According to -the proposed model, for a. flaw -speed of 70 ft/sec (21.34 
m/sec], the structural motion favorable to a compliant effect has a wavelength 
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o£ 0,066 ins. Cl *68 mQ., an amplitude* of. -660 vixns.. (16*7.6 ym)v.and a frequency 
of 3700 Hz. 

THE EXPERIMENTAL PROGRAM AT, LANGLEY- RESEARCH = CENTER 

The compliant wall experiments conducted by the Fluid Mechanics Branch, 
High Speed Aircraft Division of Langley Research Center can be. classified into 
two categories of test -conditions; (1) Wind tunnel test conducted in a Low 
Turbulence Pressure Tunnel (LTPT) with flow speeds of 210 ft/sec (64 m/sec) 
and higher. The results of. these-. experiments, have been reported (ref. 5). 

(2) Tests conducted in a 7 in.' by. JLl in. low speed ttinnel where flow speeds 
can be varied between 50 and 150 ft/sec (15.2 to 45.7 m/sec) . 

Because of the differences in test, conditions,- model, geometries are scaled 
differently for these eitperiments . So far, LTPT tests have been made in lamin- 
ated structures only and the 7 in. by 11 in. tunnel tests have been done on a 
wide variety of models. 

Currently, drag measurements are. made using a direct drag balance, and 
facilities for monitoring surface motions during an actual run. are available. 
Thus, in some. cases, comparison between the theoretical data presented in this 
report and the experimental .values has. been possible. 


OVERALL PERSPECTIVE OF IHE STRUCTURAL ANALYSIS PROGRAM 
related’ TO .COMPLIANT WALL EXPERIMENTS 


Material- Requirements 

The compliant surfaces used for testing have been built of elastomers. The 
manufacturer’s quoted properties of. the variotis- film materials as. well as the 
foam substrate are often an imreliable guide to the true conditions in the stock 
material. The compliant wall materials used in the Langley experiments were 
therefore carefully tested in the Materials Division of. Langley Research Center. 
Materials were tested in film foam, on an Eustron model TTC tensile testing 
machine. Typical compliant materials used in the experiments were PVC plastisol, 
Alathon, latex rubber, aluminized PVC, polyethylene and aluminized mylar. A 
brief summary of materal testing is reported elsewhere (see Appendix A) . 
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Structural Response Analysis 

Previous investigations of compliant wall structxiral interaction 8, 

9, 10) have been based on speculation that the compliance of the structure acts 
as a stabilizer of the transitional instability. Neither the effect of supports 
nor the presence of nonlinear effects in the motion received any attention in 
these works. In some later works (ref.: 11) the viscoelastic nature of the material 
has been -included in the analysis . However, the problem of fluid-structural 
interaction has largely been given only cursory attention. 

The present work has attempted to remedy this situation in two directions: 

1. Modelling of the actual structure, has been done carefully; wherever 
numerical simulations are appropriate they have been used. IVhen semi-analytic 
methods are useful we have employed these. Identification, of .the variotis loads 
acting on the structure tmder test. conditions is. another major. problem given 
attention in this report. 

2. Inclusion of structural nonlinearities or other, gross effects which 
can affect the actual motion has been made after proper evaluation (for 
instance, if the loading ,is moderately small, structural nonlinearity need not 
be included when the nonlinear effects are mainly geometric) . If acoustic 
effects are important these effects need to be included in the analysis. If 
mean flow effects come into the picture (flutter, divergence) then .the appro- 
priate loading due to this needs to be included. 

Because of. the wide variety of compliant models, we resorted to a classi- 
fication which covers basically most of the compliant experiments reported in 
the literature: 

1. Laminated structures: Usually membranes bonded to substrates which 

are rigidly fixed at the bottom. The membranes- may be under tension; 

2. Membranes under deep cavities l The cavities could be fluid filled; 

3. Membranes under narrow air gaps; and 

4. Periodic structures: The periodicity due to. multiplicity of equally 

spaced supports usually in the flow direction. 
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STRUCTURAL OSCILLATIONS IN . THE- PRESENCE OF FLOW 


The response of a structure to loading in the presence of. a fluid is quite 
different from its in vacuo response. The motion of the structure introduces 
a disturbance field in the fluid which itself is felt by the structure. Thus 
there are two distinct boundary value problems that need to be solved: 

1. The boundary value problem .of the structure of- volume V, with the 

given dynamic loads external to it and body forces: only a part of its total 

boundary surface S — say, S 2 — interacts with the fluid. 

2. A botindary value problem for a .fluid, region. R interacting with the 
elastic body, the solution of which has the form 

p = pCU, x) (13 

where p is the hydrodynamic pressure at the point X on and U is the 

displacement of the points of the body on the boundary. 

We thus break xjp the total loading acting on the structure into three 
components : 

a. P , the external pressure loading acting on the structure causing 
the initial or primary motion. 

b. P^, the induced pressure loading on the structure due to its inter- 
action with the ‘primary fluid. 

c. P^, the back pressure or the induced pressure from the back side of 
of the structure. 


External Force Field 

For compliant wall experiments under consideration in this report, the 
external pressure loading P can be broken down into three components: 

1. A random convected pressure .load due to the. turbulent boundary layer. 
(The pressure pulses travel at speeds between 0.5 to 0.8 times free stream 
velocity,) Bull (ref. 12) has made measurements of the spectral pattern of 
this pressure for subsonic flows. At lowest wave numbers, Bull’s spectrum 
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does not drop off, and we liave, vised the .dropoff .suggested by -Von Karmen a^d 
Lin (ref. 13}, i.e., E(k,t} == FCl,s°,k,t) where e is the eddy diffusivity, 
and I is Loitsianskii's integral. 

Wherever only qualitative trends are desired (e-g*! expected value of 
surface amplitude, power spectral. response} we have used Bull’s spectrum (with 
the proper modifications) for analysis. 

Recently, we have developed a full simulation of the random pressure field 
using Monte Carlo techniques (ref, 14} . In some anal.yses using numerical simu- 
lation for structural response (specifically membranes under thin air gaps, 
laminated structures using NASTRAN) , we have used the simulated pressure field. 
In figures 1 and 2 we display the modified wall pressure spectrum and the simu- 
lated wall pressure history at points for a low speed turbulent boundary layer. 

2. Static pressure differentials. In some wind tunnels there will be a 
difference between the static pressure of the fluid and the ambient pressure, 
which will be a function of the flow speed. A knowledge of this deviation is 
essential since it can cause primary surface deformations. 

3. Pressure gradients. The presence of a static pressure -gradient in the 
tunnel affects the theoretical predictions in two' ways. In all of the analyses 
presented here, a tacit assumption is -made that the mean flow is parallel shear 
flow. The nonparallel effect due to gradients requires a nonparallel stability 
model for the proposed model (ref. 1) . Furthermore, the induced pressure force 
evaluation also suffers from lack of consideration of the nonparallelness. 
Hence, if pressure gradients are present in a given experiment, we shall ignore 
that experiment altogether. 

Induced Pressure Field P^ 

The motion of the structure introduces a perturbation' pressure field into 
the flow and hence on itself. Sophisticated fluid theories can be developed to 
evaluate the boundary value problem, of a flow of infinite extent — one boundary 
of which changes with time. Since many of the experiments. were conducted in 
the predivergence regime we have accepted a lower order theory, viz. piston 
theory (potential flow model} , for developing the e^qiression for induced 
pressure on the structure. 


6 




P(Xg,t) (N/m^) 



The inclusion of the induced pressure field is unwarranted in some cases . 
These are cases where the divergence speed of the model is fax above the opera- 
ting speed. An analysis which neglected these forces would still yield accurate 
results for structural response in these cases. As we proceed in these analy- 
ses we shall indicate the inclusion (or otherwise) for different experiments 
justifying our reasons for same. 

1 ^ 

Back Pressure P 

Some of the models tested contained fluid filled cavities beneath them. We 
shall develop theoretical expressions for evaluating^ these in the next sections. 

ANALYSIS OF A MEMBRANE OVER A DEEP CAVITY WITH EXTERNAL 
FLOW OVER THE MEMBRANE 

We consider (fig. 3) a membrane of flexural rigidity. D, tension T, 
simply supported on ends x = 0, x = a, and y = 0 and y = h, and of thickness 
h and density 



Figure 3. Membrane over a deep cavity. 
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The dynamic equation of motion for this membrane is given by 


3t^ 


C 1^ ? DV^V^w 
at 


- TV^w = P 


e . 

+ P 


+ P 


C2) 


In equation (2) we have neglected effects of midplane stretching. IVhen non- 
linear Cs'txetching} effects need to be taken into consideration, we use the 
following equations: 


ph ^ + C IJ- + DVV - ( 

9t2 . V 


N ^ + 

ax^ 


N 


a^ 


w 


+ .N 


a^w \ ^ 


y 3 „ 2 . ' xy axay/ 


p® + + p^ 


C3) 


and 



ph 




1 - v\ a^v 1 + V a^u 
, 2-/ 3^2 2 8x3y 


(4) 


a^tf / a^w 1 - V 3^w\ 1 + V a^w aw | _ 

3y\ay2 2 g^ 2 j 2 axayaxj=" 


and 


N 


x 


N 


y 





cs;) 


C6a) 


(6b) 




3v ^ 3w 3w \ 
3x ax 3y I 


(6c) 
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Linear Model 


We shall assume that the surface motion is made of normal modes in space 
and has the form 


CO CO 

w = V W (t) sin 
^ ^ mn 
m n 


— sin 
a b 


( 7 ) 


Furthermore, we assume that a harmonic dependence on time exists, i'.e.. 


W (t) = a e 
mn ^ mn 


lU)t' 


(8) 


For the linear problem at hand, equation (8) is a justifiable assumption 
as long as P® can be broken down into its harmonic components, i.e.. 


ioj.t 


p" pj = ' 


Evaluation of the Back Pressure 


Under the assumption that the fluid motion in the cavity is a perturbed 
state of motion and is irrotational , we define expansions for the velocity 
potential and the elevation of the surface Sx as 


ij, = £ + e2 ^(2) ^ _ 

, = , ... 


03 




(X) 


satisfies the wave equation 


, 2 ^, 0 =) 

c2 3t2 


CIO) 


where c is the speed of sound in the medium. 


11 



The boundary conditions are, at rij i.e., on 

[t) = z - w(x,y,t)] 

= n - ^ 3'f’ 

dt ~ ° ~ 3t 3x 3x 3y Sy ■ 3z 


and at z = H 


91 

3n 


0 , 


i.e 


♦ i 


3(j) 

3z 


0 


Applying the expansion, equation (9] into equation (11) 


3z 


( 0 ) 


at 


= 0 


3z^l^ _ 


at 


3z 


yielding 


CO) 


= constant = 0 , 


and since 




= 11 I 

3t 3z I 


z=0 


( 11 ) 


( 12 ) 


(13) 


Let us now assume that the solution of equation (10) with the boundary condi- 
tions [eqs. (11) and (12)] has the form 


<{. = a(z) a^^(t) sin-^,sin-^ 


(14) 
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Thus 


a"Cj) - 1(E) + } atz) = ^ aCz) 

with a' CO) = iti) 
a’C-H) = 0 


Putting 



(IS) 


the solution for is given by 




io] 


cosh 


X 


mn 


X (z - H) 

_ mn^ ^ 

sinh(X H) 

^ mn 


a (t) 


sin 


mirx . mry 

sin 

a b 


(16) 


In order to evaluate the back pressure, P^, we apply Bernoulli’s equation 
at z = 0, 

gw. 11=^ C17) 

b 

or 

0 )^ coth X H 
ran 

X 

mn 



. mTTX . mry 

a (t) sin sin -~ 

mn^ a b 


(18) 


Evaluation of the Induced 
For the velocity potential we 

c^ ' 


Pressure Due to the Primary Flow 
can write the governing equation as 


(19) 
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From momentum balance in z direction for the perturbed Eulerian flow 


p = -p[^ * 


(203 


The boundary conditions on (() are 


9i() I 

3z I 


= 0 , 


3<l>l / 3 . .r 9 ^ 

z=0 ~ ’ 


( 21 ) 


using equation (20) we rewrite equations (19) and (21) as 




( 22 ) 


with 


3P 
3z ' z 


-0 ■ ” ’ 


C23) 


3^1 

3z ' z=<» 


== 0 . 


(24) 


Solution for 

Once again we use for w the form given by equations (7) and (8) . Because 
of the nature of the boundary conditions [eqs. (23) and (24)] and the governing 
equation (22) itself, we can seek the solution p at the wall z = 0 for a 
mode (call it Pjjjj^^) ^.nd using superposition obtain 

P (25) 

mn 

m n 


We shall therefore consider the expression for w as 


, , not . mvx . mry 

w = a e sin — — sin -r^ 
mn a b 


(26) 
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and rewrite 


w 


a 

mn 

4 


[Wx + W2 + W3 + W4] 


and we i\rrite 


w 


[wj + W2 + W3 + W4] 


C27) 


(28) 


where 


w. 


i(k)[x cos 0. + y sin 0. + Cit] 
e 3 3 


j = 1,2, 3, 4 


(29) 


and 


e.'s 


are such that 


01 = -0 

02 = IT - 6 

03 = 0 

04 = TT + 0 


(30) 


where 



(31) 


k 



(32) 



(33) 


and 


Cl 


0} 


oja 


mfrg 


mn 


(34) 
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In order to evaluate the induced pressure at the wall we compute the pressure 

due to w. . 

1 

The solution of equation (2) with equations (23) and (24) for an input 
wall function w is given by 


P(z,x,y,t) a® 


mn 


-|klp^z 


m 


n j=l 


i|k| [x cos e. + y sin 9. - cjt] 
a I 1 


(35) 


where 




'U cos 6 . + Cl 
1 _ I 1 




a . = 


V /Ci + U cos 6.\^ 
,.( 


mn 


* • — ^ 
1 4 


(36) 


(37) 


with o-i 1, a£ ~ Ij tt 3 - -1, and U 4 — -1. The wall induced pressure iSj 
therefore, given by 


“SEE 

n j=l 


ra 


ijkl[x cos 0. + y sin 8. + cit] 
a . e ^ ^ 

J 


(38) 


Evaluation of the Divergence Speeds for Subsonic 
Flows IVhere M « 1 

We shall investigate the nature of the induced pressure for sub- 

sonic flows with M << 1: 

1. For 03 = 0 (for m = 1, n = 1) we obtain 


i 

P 


a X 
mn 


Ikl^ 


A ~ 


sin 


m-irx 

a 


sin 


niry 

b 


(39) 
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2. For nonzero w 


p = P 1 + i P7 

real imaginary 

The imaginary term acts as a severe damping term. 

The critical speed can therefore be defined as the speed at which the 
membrane goes into zero frequency flutter. Since frequency is zero, the 
lowest critical speed for the structure occurs when the induced pressure p^ 
is equal to the stiffness of the structure. Thus, 



(40) 


For a given ratio of I?; the flutter boundary (i.e., critical speed) can 
be found from equation (40) . 


Nonlinear Effects 

The linear model is only capable of specifying the divergence speed. In 
an actual experiment one would see membrane divergence and small frequency 
flutter around the critical speed. This is because the nonlinear effects 
have not been included in developing equation (38); also, the critical speed 
will be slightly higher than that predicted by equation (38) since the effect 
of nonlinearity is to increase the stiffness. In order to obtain a valid 
analysis one has to: 

1. Assume that a^^ are slowly varying functions of time, and conduct 
an analysis using Kelvin's stationary phase method to obtain closed form 
solutions; or 

2. Develop numerical analytic solutions from p based on Fourier 
transform techniques (see ref. 15) . 
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STRUCTURAL ANALYSIS FOR LAMINATED STRUCTURES 


In many compliant experiments reported in the literature (refs, 4 and 5)^ 
the compliant model consisted of a very thin membrane (thickness ~ 1 mil) 
bonded to a soft foam (E ~ 200 psi) of thickness 0.25 to 0.5 inches. The 
analysis of the "laminated" structure is performed using proper structural 
modelling. In the following we present different approaches in modelling 
the structure and their basic merits. 


Approximate Analytical Model 

Consider (fig. 4) a membrane (a x b x h) with properties (P^, E^, 
under uniform tension T which is bonded to a substrate (a x b x H) of prop- 
erties (P^, E , V ). We assume that the substrate can be modelled as a semi- 
s s s 

infinite foundation undergoing plain strain deformations to external loading 
p(x,y,t). The end effects near the edges are thus neglected. The effective 
foundation properties are 



Of the load pCx-,y,t) acting on the membrane we assume that a portion q(x,y,t) 
is transmitted to the foundation. The deformation state of the foundation is 
therefore assumed as 


u(x,y,z,t) = v(x,y,z,t) = 0.0 
w(x,y,z,t) = W(x,y,t)ij)Cz) 


(42) 


where W(x,y,t) is the deformation of the midplane of the membrane. 
The function ^^(z) is such that 


ii>(0) = 1.0 

ijj(H) = 0.0 


(43) 
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Figure 4. Laminated structure. 





The functional form chosen for ij; is 


ilj(z) 


sinh u(H - y) 
sinh uH 


( 44 ) 


where u is a foundation constant. The stress components -of the foundation 
are 



3w 


4' 


8W 


'zy 2(1-5- vp dy 2(1 + ^ dy 


dw 


zx 2(1 + v^) 3x 2(1 -I- v^) 


3x 


(45) 


Applying the principle of .virtual work^ the equation of motion of the substrate 
under the given loading q(x,y,t) is obtained as 


/ 


H 




ifidz 



a III 'dz 
z^ 



— ijj^dz = -q(x,y,t) 
at^ 


(46) 


Equation (46) can be \tfritten as 


q(x,y,t) 


m* 


k*W - T*v2w 
at2 


(47) 
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where 


m* = p 


-H 

''0 


Cl - v|) *^0 


/ [r]^ dz 

r\ 


: /*H 

+ V J J ^ 


2(1 + V.) 


*'0 


The equation of motion of the membrane can now be developed as 


P h - — + DV^^W - TV^W = p(x,y,t) - q(x,/,t) 
® 3t2 


where 


flexural rigidity = 


12 Cl - V-) 


or, using equation (48), 


(p h + m*) + DV*^W - (T + T*)v 2W + k*W = P(x,y,t) 


The eigenvalues of the structure can be determined from equation (51) as 
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One o£ the chief merits of an analysis such as the one above is its simplicity. 
However, the fidelity of the solution so developed depends strongly on the 
correct choice of the foundation parameter u. 

Models Based on Simulation Techniques 

An implicit assumption in the preceding section was the one-dimensional 
(depthwise) variation of all deformations. The foundation thus acts basically 
as a series of vertical springs supporting the membrane. The finite nature 
of the structure makes the approximation carried out in that section rather 
inaccurate. A proper approach to modelling would be using k numerical 
technique, viz. finite elements. 

Using NASTRAN (acronym for NASA Stru ctural. Analysis Program) we per- 
formed simulation studies on the laminated structure. The basic models 
are briefly described below. 

Two-dimensional model . The structure is assumed to be two-dimensional 
(X-Z) in the direction of the flow (fig. 5a). The membrane is modelled as 
bar elements (CBAR) with only flexural properties. The substrate is assumed 
to develop only shear deformations due to the external loading] hence, quad- 
rilateral membrane elements (CQDMEM) are used to model that. OFFSET cards 
are used to account for the offset of a node of the membrane element of the 
substrate from the node of the bar element at the interface. 

Plate-spring model . A static three-dimensional analysis is made for an 
elastic slab (substrate) using three-dimensional isoparametric elements 
(CIHEX3) (see fig. 5b). From the static analysis the equivalent spring 
stiffnesses of the foundation are obtained. The membrane is modelled as 
a plate with membrane action (CQUADl) . For a given loading the structure 
(CQUADl-CELAS) is analyzed using NASTRAN. The inertia effects of the sub- 
strate are lumped to grid points using C0NM2 cards. 

Fully three-dimensional model . The substrate is made up of three- 
dimensional isoparametric elements (NASTRAN level 15-9) CHEXA2, and the 
membrane is modelled by CQUADl. The offset between the plate grid points 
and the surface grid point of the three-dimensional elements is neglected. 
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MEMBRANES UNDER NARROW AIR GAPS 


One of the major considerations in developing compliant models for tests 
in the new 7 x 11 wind tunnel has been the ability of the surfaces to undergo 
small wavelength motions. Such motions are "favorable” for drag reduction 
according to the criterion of the Bushnell mechanism (ref. 1) . In order to 
facilitate small wavelength motions, a membrane under narrow a narrow air gap 
was tested at the Langley 7 x ll tunnel. A drag reduction of about 10 per- 
cent was obtained in this test. The compliant model is shown in figure 6a. 

In figure 6b we indicate the development of small wavelengths due to chopping 
at the bottom. In the absence of the narrow air gap the deflected profile 
would have been as indicated in figure 6c. 

The governing equation of motion for a membrane under a narrow air gap is 


ph 


+ DV^V^w 

8t^ 


TV^W a PCx,y,t) 


with 


w(x,y,tj a -6 I 


3w 

3t 


(x,y,t) 



ti _< t < t2 


where 




< 0 


and 


C53) 


CS4) 


t2 - + e 

where e is a very short duration of time whose value can be obtained 
through analysis of a dynamic Hertz contact problem. 

The motion of the membrane' under these conditions is highly nonlinear. 
Hence, no close form solutions are possible. A finite difference solution 
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for equation (S3) is possible, and an algorithm with leapfrog differencing 
.in time and center differencing in space was used to solve equation (53). 
The C-F-L limitation imposes 


Ax 

At 


> C 


sh 


( 55 ) • 


For analysis, the bending rigidity of the membrane is neglected. An appro- 
priate choice for At has been made as 


At = 


Ax 


5C 


sh 


where 


.C . = 


'sh 


ph ; 


(56) 


For analysis Ax was chosen equal to Ay (uniform mesh spacing), and the 
value of Ax was chosen to obtain resolution of the smallest wavelength 
possible. The numerical model consisted of a 10 in. x 10 in. membrane with 
an air gap thickness, 6 = 10"‘*’ in. The smallest wavelength we wanted to 
capture was of length X = 2 in. Thus, the mesh spacing was set at Ax = 0.2 in. 

A one-dimensional analysis of the same problem is possible with greater 
resolution of wavelengths. Here one would assume variations in the y direc- 
tion to be unimportant. However, the one-dimensional analysis can only show 
trends, since the nonlinear nature of ' the problem prohibits comparisons 
between the one- dimensional analogue of equations (53) and (54) , and the 
two-dimensional problem. 

The conclusions to be dra^ra from extensive simulation studies conducted 
are that (1) the model equations (53) and- (54) have unique solutions (con- 
vergence with reduction of step size) and (2) the solution for dynamic 
response exhibits small wavelength, high frequency motions compared to a 
pure membrane. 
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PERIODICALLY SUPPORTED STRUCTURES 


Compliant drag reduction has been reported in the literature for period- 
ically supported structxires (ref. 2) . The length of the bays for these com- 
pliant surfaces corresponds to the wavelengths predicted by the Bushnell 
mechanism. We develop below an analytic method of solution for periodic 
structures based on the work of Mead and Pujara (ref. 15), 

The theory developed here is basically a linear small deflection theory 

with the requirement from consideration that the midplane strains are negli- 
w 

gible, i.e., << !• We will therefore discard solutions which give ampli- 

tudes w ~ 0(h) or greater. The above requirement is essential, as the ratio 
of the linear restoring force versus the nonlinear stretching force needs 
to be very large, i.e., 


w X 


Eh 3 


12(1 




(¥) 


- » „3 Eh(f) 


or 

1 

h2 


(57) 


We also assume a priori that the structure is one-dimensional, i.e., crossflow 
directional variations are neglected (a/b « 1) * 

Figure 7a shows a periodically supported beam, over which is a flow with 
speed U^. Figure 7b shows an equivalent representation for flexibility at 
the supports and sitting on an elastic foundation (stiffnes = k^) . If the 
foundation is viscoelastic the representation for the foundation is made 
= k^(l + i n^) where is a typical "loss factor" for the foundation 

and k^ is the complex stiffness of the foundation. 

We shall develop the theory of response of the structure for two cases: 

1. The foundation is absent; instead there is a cavity of fluid at the 
backing, and 

2. The case with the foundation. 


27 








Case 1. To develop expressions for this case one has to evaluate the 
reaction of the cavity on the bay under consideration. The approach 
is very similar to evaluation of the induced pressure acting on the 

side where there is flow and hence will not be dwelled upon in detail here. 
Also, one has to set the foundation reaction (k^ = 0) to zero. 


Case 2. Consider an external force field 


-i!i ^ 

P = P e e 

0 


CS8) 


We seek harmonic solutions to the excitation. The transverse displacement 
can be written as 


w(x) 



A 

n 


-iCp+2nir)j 

e 


i(i)t 

e 


(59) 


We shall assume that the flexural rigidity of the beam is D’ 


D’ = D(1 + in^j) 


(60) 


Eh^ 


where n is the loss factor for the beam and D = - 

12(1 - v2) 

of motion of the beam under consideration of excitation is 


The equation 


D' 


— + p,h — + k' w = P®(x,t) + P^(x,t) + Ri 6Cx) -i- R 2 6Cx - £) 
Sx*^ ° 3t^ ^ 


( 61 ) 


where R^ and R 2 are the reactions of the supports and P is the induced 
pressure. The induced presstire can be evaluated as was done in previous 
sections as 


p^(x,t) 
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(62) 


where is the density of the fluid medium, 


29 




and 'a’ equals the speed o£ sound in the medium, 

+ 2nir (64) 


In order to obtain solution of equation (61) for given loading [equation 

- i(ym - mt) 

(.58)], we evaluate the virtual work done by virtual displacement e 

and equate the sum of all virtual work to zero. Because of the periodicity 
of the bays the total virtual work done by the supports to be included in a 
single bay is just due to that of one of the supports. 


The virtual work done by the support at x = 0 is 



6Am 1 K ' Y A + K 
t ^ n r 
^ n=-«> 




(65) 


The virtual work done by the other reactive forces can be obtained by multi- 
plying equation (61) by the virtual displacement and integrating over the 
interval 0 to 5,. Finally, we indicate the set of simultaneous equations 
obtained by setting the sum of all virtual work to zero as 
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= 0 when m 5^ 0 


= p when m = 0 


( 66 ) 
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The simultaneous set of equation (66) can be solved finally to obtain 
A., In practice, only a few of these A^'s have to be evaluated. Taking 
n to be between -6 to +6 is sufficient to obtain the desired accuracy. 
Usually the rotational rigidity of the supports is zero (simply supported, 
i.e., Kr = 0) we indicate here the approach which is simple to solve 
equation (66) for this case. 

We rewrite equation (66) as 


A„<{>+4^Za=P m?^0 

m^m Z ^ n o 

= 0 m = 0 


(67) 


where 
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From the case m = 0 Ave obtain 
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for the case where Kp equation (69) yields in the limit 


A = 
m 



CO 




-CO 



(70) 
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The evaluation of leads to the knowledge of w(x) for a given harmonic 

excitation. 

In order to obtain the solution for a wide band excitation one must carry 
out a Fourier analysis of the exciting field and stun up, using the principle 
of superposition, the harmonic responses given by equation (59]) . The wide 
band response can then be analysed to obtain power spectral data of expected 
response. The mean square response is defined as 

E[u)0)*] = Ejp(uj)P*(a))} (71) 

where H(j&}) is the frequency response function and e|p(oj)P* (( u)| is the 
power spectrum of the turbulent pressure field. 

NUMERICAL RESULTS 

In the preceding sections we have developed appropriate techniques for 
analysis of a generic of compliant structures. Using these analyses we 
examine in this section the nature of structural response under reported 
experimental conditions . 


Laminated Structures 

The compliant model tested in the LTPT tunnel at Langley consisted of 
an 51 in. x 18 in. model surface. The compliant surface was made of 1 rail 
mylar (Young's modulus ~ 5 x lo^ psi, v = 0.3) with a backing substrate 
which is a soft foam (80 PPI Polyurethane foam, E = 200 psi) of thickness 
0.25 in. In table 1 we report the eigenvalues (natural frequency of vibra- 
tion) obtained using different models discussed in the Structural Analysis 
for Laminated Structures section. These eigenvalues are arranged in ascend- 
ing order. Those eigenvalues 'listed in the last column under "Membrane" 
represent the natural frequency of the membrane (in vacuo) if the substrate 
\-jeve absent. Figure 8 presents the frequency response of the laminated 
structure and indicates that the amplitudes are small; the turbulent pressure 
excites the structure into a wide band response with most of the power 
included in the frequency range below its lowest natural frequency. 
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Table 1. Natural frequencies of a laminated structure 


Eigenvalues 

Using NASTRAN 
Models 

Using Eq. ’ (52) 

Membrane 
In Vacuo 

C1,1) 

484.6 Hz 

535.2 

125.6 

C2,l) 

485.0 

536.0 

145.0 

(3,1) 

485.7 

537.6 

172.5 

(4,1) • 

486.6 

539.5 

204.9 

(5,1) 

486.6 

541.8 

240.2 

(1,2) 

487.7 

541.8 

240.5 

(2,2) 

488.9 

542.6 

251.2 

(3,2) 

490.4 

543.9 

268.0 

(4,2) 

490.5 

544.6 

277.4 

(6,1) 

492.0 

544.9 

290.0 


The membrane is 1 mil mylar, substrate 0.25 in. thick, 80 PPI PU foam 
(E = 2.2 psi, V = 0.1). The model is 51 in. x 18 in. and with a tension, 
on the membrane of 2 Ibs/in. 


** The foundation parameter v) is chosen to be 2.0 















We conclude that laminated structures do not show any promise as "passive" 
compliant surfaces on the basis of the Bushnell mechanism.- 

Membranes Under Narrow Air Gap 

As was indicated in the Membranes Under Narrow Air Gaps section above, 
the analysis of the membrane over a narrow air gap only presents qualitative 
information regarding the nature of the motion. Using the simulation program 
for pressure (ref. 14) we have developed solutions for a 10 in. x 10 in. 
compliant model with 1 mil (1000 pin.) mylar over 100 pin. air gap and a 
flow of 50 ft/sec. Such flow conditions typically occur in a 7 in. x 11 in. 
tunnel at Langley. The boundary layer thickness was approximately 0.5 in. 

The pressure history and displacement for the midpoint of the surface is 
indicated in figure 9 for a duration of time. 

Figure 10 indicates the influence of spatial discretization on the solu- 
tion obtained. Figure 11 presents an analysis of the frequency response of 
the motion. It is shown that the narrow air gap-membrane configuration 
yields high frequency motion. 

Membranes Under Large Cavities 

In some reported experiments (ref. 3) membranes over fluid-filled cavities 
have shown drag reduction. In order to understand the nature of the surface 
motion we examined a recent compliant wall experiment: the membrane is 1 mil 

thick mylar and the cavity is 0.25 in. deep and filled with air. The test 
speed is 50 ft/sec. The mean square response with frequency is indicated 
in figure 12. It is shown that the motion is purely low frequency, large 
wavelength. Therefore, we do not believe that membranes over large air gaps 
are promising candidates for compliant experiments. 

Figure 13 indicates the response of a membrane with a water-filled 
cavity. Again the presence of the cavity shifts the response curve towards 
the low frequency end of the spectrum. 
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Figure 9. , The pressure history and resulting surface motion 
for a membrane over an air gap (U^ = 15.2 m/sec, 
2.54 cm boundary layer thickness). 
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Figure 10. Influence of resolution on predicted surface motion. 
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Figure 11. Membrane over narrow air gap- -frequency response curve. 






Periodic Structures 


The original reporting of compliant drag reduction (ref. 2) was for a 
periodic structure made of a rubber diaphragm (80 mil thick) supported on 
periodic stubs of rubber with water filling between the stubs. The height 
of the stubs was 40 mil and the width 20 mil. Using the above data we 
examined the surface response to turbulent flow (medium is water) at a speed 
of 50 ft/sec corresponding to the experiment of Kramer. Figure 14 indicates 
the surface response. We determined that the periodic Kramer surface was a 
high frequency passive surface and could have been a drag-reducing surface 
satisfying the conditions of Bushnell criterion. In nondimens ional units 
our analysis indicates for the Kramer surface a X'*’ = 800, h"*" = 100, and 
c = 0.7 U^. The analysis includes a nominal value of damping for the 
structure. We expect, however, that the actual damping in the structure 
may limit the amplitude to h — 20 or 30 in the experiment. ' Such an exci' 
tation is capable of reducing the rate of burst production and hence could 
have been responsible for the observed drag reduction. 

Most attempts toward reproducing the results of Kramer over the years 
have been unsuccessful. One of the chief reasons for this has been the 
belief that Kramer's experiment delays transition; hence, subsequent experi- 
ments at various laboratories were conducted as transition experiments. In 
a few cases the flow had been accidentally tripped to turbulence, and there 
have been reports that drag reduction was indeed observed for these cases. 
However, the majority of data points showed no favorable drag changes (ref. 
17) . We also emphasize that proper attention was not given to the structural 
motions required. Consequently, the geometries tested were not facsimiles 
or scale models of the Kramer surface. There is, therefore, a need for re- 
examining Kramer's surface for possible drag reduction. 

PASSIVE WALLS FOR 7 IN.^x n IN. TUNNEL 


The 7 in. x n in. tunnel at Langley has a speed range of 50 to 150 
ft/sec. The flow is tripped about 1 ft from the test model using rougheners 
(sandpaper) . Boundary layer surveys on the test model indicate that a rela- 
tively thick boundary layer (6 = 0.5 to 1.0 in.) is formed on the model. 
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Figure 14. Steady- state response of Kramer stubbed periodic surface 
= 18 m/sec water flow) . 



The model itself is 11 in. x 36 in. and floats on a drag balance (see fig. 15). 
The facing plates smoothly fan the flow onto the test area. IVhen rigid plate 
measurements are made the whole surface area is available. However, the 
compliant section on which measurements are made spans only a 7 in. x 16 in. 
area in order to avoid nonuniformities of flow around the side walls. The 
adjustable backpiate provides a facility for adjusting the cavity depth which 
can be filled by a foam or with any fluid. IVhen a backing such as air is 
used, the back pressure is controlled using a vacuum pump. The overall 
pressure gradient in the test region is kept at a minimum. However, there 
is a slight pressure gradient of a few y psi/in. in the tunnel test section. 
The total pressure on the flow side fluctuates slightly from the ambient, 
and the tunnel static pressure differs from the- ambient pressure as a small 
function of the flow speed under conditions of optimal operation. IVhile 
evaluating the response of structures under flow we included this small 
static pressure differential acting on the structure in addition to the 
existing force fields. 

The vacuum section chamber in the model (fig. 15) is used to put tension 
on the compliant wall. The widths of the cavity of the suction chamber are 
adjusted in such a way that the tension is uniform. This gives the ratio 
of the cavity widths in x and y directions as 



(see Appendix B) in order to obtain uniform tension of the membrane. By 
directly observing the center deflection of the membrane in the cavity on 
any side and noticing the pressure differential, one can directly evaluate 
the tension in the membrane knowing the properties of the membrane. 

Finally, direct calibration can be made between the tension in the mem- 
brane and the vacuum pressure provided that the properties of the compliant 
membrane are fully known. 

In figure 16 we show the observed surface motion using a schlieren- system 
for taking area photographs. The surface is a membrane over a narrow air 
gap. The dark spots in the picture correspond to points where the membrane 
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Figure 15. The compliant wall model. 
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Figure 16. Area photograph of the response of a membrane 
over a narrow air gap; backing surface is 
10 PPI foam. 
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touches the 10 PPI foam surface after the excursion through the narrow air 
gap. The arrows indicate the direction of flow. The tension in the membrane 
is low. 

In figure 17 we show the surface motion as a function of time at a point 
on the surface as seen in the oscilloscope. Each division in the screen 
corresponds to 5 milliseconds (in the x direction] and in the ordinate 
1 division corresponds to 0.1 volts of photodetector output for the bottom 
picture and 0.01 volts for the top one. The output of the photodetector is 
directly proportional to the maximum angles of the surface motion and has 
been calibrated already. 

From figure 17 we learn that the surface motion of the membrane over the 
narrow air gap is high frequency (between 400 Hz and 2 KHz) and has ampli- 
tudes of 100 pin. or more. 

In figure 18 we show the area photograph of a "dulcimer” model. The 
model consists of a membrane sitting on periodically supporting strings. There 
is a separation between the membrane and the strings, and during surface 
oscillation the membrane slaps on the strings. The dark lines correspond 
to the strings in contact with the membrane. The strings are kept in tension. 
The motion of this surface is again small wavelength and at moderate fre- 
quencies. Figure 19 shows the motion of a point with respect to time. The 
time scale is 20 milliseconds/division. The“flow speed is U = 160 ft/sec. 
Again, the surface amplitudes are only about 200 pins. (1 voltsslOOp in.). 

Using the theoretical analyses in the preceding section-, we also con- 
ducted some experiments where the wall was made to flutter. (The spacings 
between supports were designed to be such that low' frequency flutter occurred.) 
One of the chief problems in these experiments was adjusting the pressure in 
the back chamber continually while the experiment was going on, in order to 
suppress the panel divergence. 

The models were tested for drag reduction too. In most of these compli- 
ant experiments there was little drag change. The observed motions did not 
suggest that the surface motion was in the range of parameters suggested by 
the Bushnell mechanism either. Further compliant experiments where compli- 
ant motions suggested by the Bushnell mechanism can be produced (active walls/ 
passive walls) are unden\?ay. 
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Figure 17. Surface motion of a point on the membrane over 
narrow air gap as function of time. 
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Figure 18. 


Area photograph of the surface motion of a 
"dulcimer" model . 
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Figure 19. Surface response with time for the dulcimer model. 


48 



CONCLUSIONS 


Theoretical analyses for common compliant wall geometries are presented. 
Analyses of surface motion for the compliant experiments indicate that the 
original Kramer experiments conducted in water could have' interrupted the 
turbulent burst mechanism and produced the drag reduction. Membranes over 
thin air gaps can also be used to produce the small wavelength, large fre- 
quency motion required for air experiments. It is possible that in some 
laminated surface compliant wall experiments, delamination of the structure 
occurred and thereby created a narrow air gap situation; thus the observed 
drag reduction in the experiment of Walters (ref. 4) and the Langley experi- 
ment (ref. 5) . The membranes over large air gaps do not sho\tf any promise 
at all. The "observed" drag reduction reported in reference 3 can be due to 
other effects, as suggested in reference 16. It seems rather difficult to 
obtain using passive walls the desired surface motion for most low speed and 
transonic speed experiments in air. Thus it may well be worthwhile to use 
active walls, i.e., walls where the motion is driven onto the surface for’ 
compliant experiments. Finally, there is a clear need to understand the 
entire fabric of fluid structural interaction, and this can only be achieved 
by bringing in and coalescing structural expertise with the fluid mechanics 
of turbulence. We believe that this report goes as far as making a first 
tentative step in the right direction toward such a goal. 
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APPENDIX A 


FILM MATERIAL TESTS 


By 

Todd Hodges^ 


INTRODUCTION 

Tests have been completed on eight materials in film form to support the 
compliant skin drag reduction project in the Fluid Mechanics Branch. The 
Fluid Mechanics Branch was interested in the elastic modulus of the materials 
they were using in their wind tunnel models. 

TEST EQUIPMENT AND MATERIALS 

Tests were conducted on the Eustron model TT-C located in room 003 of 
building 1293A. The C range load cell was used. Depending on the material 
being tested, the full scale ranges used were 1 lb, 2 lb, and 5 lb. Accuracy 
of the load cell is equal to, or exceeds, 0.25 percent of the range in use or 
0.50 percent of the indicated load, whichev er is greater. The amplifier and 
recorder accuracy is within 0.50 percent of full scale range. 

Eustron air grips were used to hold the test specimens , Test gage length 
was 3 inches, and specimen width was 1/2 inch. Film thickness was measured 
by a power-driven film micrometer (Model 549E by Testing Machines, Inc.]. 

The test specimens were mounted and pulled at a constant crosshead speed 
of 1/2 inch per minute. The tensile modulus (E - psi) was calculated accord- 
ing to the attached equation derivation with an electronic calculator. Five 
samples were used on each test run except for the Latex rubber, where only 
three samples survived the film cutter. 


^ Research Associate, Virginia Polytechnic Institute and State University, 
Blacksburg, VA. 
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DISCUSSION 


The following modulus values were measured on the eight systems’: 

1. Ethylene vinyl acetate copolymer (Alathon) (tested in the longitud- 
inal direction). E = 20.0 x 10^ psi with a high of 20.6 x lo^ psi and 

a low of 19.8 X 10^ psi. This material showed good consistency in thickness 
and modulus . 

2. Alathon (tested in the transverse direction). E = 15.5 x lo^ 

^ average 

psi with a high of 16.5 x lo^ psi and a low of 13.1 x lo^ psi, also showed 
good consistency in thickness but not as good in modulus. 

t 

3. Latex rubber, ^average ~ high 313 psi, low 272 psi. 

Lowest modulus of all the materials . We were able to salvage three test 
samples from the film-cutting device. 

4. Aluminized kapton, Eg^ygj.age ~ ^ high 468. x lo^ psi, low 

390 X 10^ psi. Modulus consistency was not good. 

5. Aluminized P.B.C, E = 361 x 10^ psi, high 416 x 10^ psi, low 

296 X 10^ psi. Thickness and modulus variation were high. 

6. P.V.C.j E = 33.4 X 10^ psi, high 36.9 x IQ^ psi, low 28,2 x lo^ 

psi. Thickness and modulus variation were high. 

7. Polyethylene, = 2.17 x 10^ psi, high 22.8 x lo^ psi, low 

20.8 X 10^ psi. Thickness and modulus consistency were good.’ 

8. Aluminized mylar, E = 451 x iq 3 psi, high 536 x lo^ psi, low 

SCV 

367 X 10^ psi. Thickness and modulus consistency were not good. 

The aluminized films generally showed less consistency in thickness and 
modulus than plain films. This could be due to the coating technique they 
use and its effect on the substrate film. 

Ranking the materials in order of modulus we have: 

1. Latex rubber: E = 296 psi. 

av ^ 

2. Alathon (transverse): E = 15.5 x lo^ psi 

3. Alathon (longitudinal): E^.^ = 20.0 x 10^ psi 

4. Polyethylene: = 21.7 x lo^ psi 
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5. P.V.C.: = 33.4 x 10^ psi 

6. Aluminized P.V.C.: E = 361 x lo^ psi 

7. Alumnized kapton: E^^ = 390 x iq 3 psi 

8. Aluminized mylar: E = 451 x lo^ psi. 

3.V * 
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APPENDIX B 


VACUUM SUCTION TECHNIQUE FOR TENSIONING COMPLIANT MEMBRANES 


In this section we shall briefly discuss the vacuum suction technique 

for tensioning the compliant membranes. The cross section of the vacuim 

suction chamber is shown in figure B-1. The shaded area of (L ^ L ) is the 

X y 



Figure B1 . 


working area of the model on which a uniform tension needs to be applied. In 
order to apply uniform tension in the working model the width of the cavities 
£ £y have to be related in some fashion to the dimensions L , L 

Initially the membrane is taped to the sides ABCD with care such that 
there is little tension on the membrane. The suction chamber pressure is at 
ambient pressure during this time. Now the chamber pressure is adjusted to 
be lower than the ambient pressure by a magnitude Ap. The deflections of 
the membrane at the center of the channels wi, W 2 are noted. 


In figure B-2 we indicate schematically what happens to a cross section 
of the membrane. For the simply supported membrane of length £ and width 

A 2C 

(L^ + 2£y) where the ratio 
equation as : 


X. \ 

^ , I « 1 we can approximate the governing 
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one obtains for after substituting in equation (B-1) 
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Thus ail odd modes are present in the solution and 
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Similarly 


W2 


0.125 



Ap 


Thus 


CB-6) 



Now, the stretching that the membrane undergoes in the cavity 2 is con- 
tributed by the stretching o£ the membrane in the working area. The stretch- 
ing in the cavity 2. is given by 



dx - 2 
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r 


m^ 


2 

0.0417Ap 

'^x 


(B-8) 


The stretching A2 is assumed to be caused by a uniform loading P per unit 
width applied to the membrane of the working' area of the model, in which case 


PL 

X 

Eh 


= A2 = 0.0417AP 
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or 


P = 




CB-9.) 


Similarly the loading on the end in the y direction is 


Z 

y 

L 

y 


0.0417 


EhAp 
T ^ 

y 


(B-10) 


Finally, we would expect to have P = and T = T for uniform tensioning 

X y 

when 



CB-ll) 


vfhich is as 


given in equation (72) . 
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APPENDIX C 

CODE FOR ANALYZING MEMBRANE WITH BACKING, 
OR WITH CAVITIES BEHIND IT 


57 



PROGRAM MA I N ( I NPUT . OUTPUT » TAPE5= 1 NPUT ♦ T APE6=0UTPUT ) 
complex DU9X2»X3 

DIMENSION AK <75*24 ) ♦ AFACT ( 75 ) *BFACT (24 ) 

dimension vn cssoo ) 

DIMENSION OMC (75,24 ) tCMl (75* 24 ) 

DIMENSION FF(aO) 

c 

c 

c 

c this program EVALUATES THE RESPONSE OF A MEMBRANE WITH DEEP , 

1 CAVITY(AIR)* A MEMBRANE WITH WATER BACKING, AND AMEMBRANE 

2 ON elastic foundation 

1 CAVITY (AIR)* AMEBRANE 

C BULL-S SPECTRUM IS USED TO OBTAIN THE EXCITATION FIELD 

CX-* INDUCED AERODYNAMIC EFFECTS USING PISTON THEORY. 

C** ■ effects of backing using POTENTIAL FLOW THEORY. 

Cit* THEORETICAL MODEL USES LINEAR THEORY 

PROPER VALUES OF DAMPING INCLUDED USING LAB ESTIMATES 
C OUTPUT — -frequency* MEAN SQUARE FORCING FUNCTION *MEAM SO 

JUARE RESPONSE FUNCTION* MEANSQARE AMPLITUDE* 

C ALL DATA TO BE IN IN-LB-SEC UNITS 




c 

c 

c 

NR = 5SNW = 6$G=386.0:eePS = 0*00000000 1 $00= 1 . 35* 1 0 . *-«-4SP 1=3*141 5926 
00 110 J=1 *5 

READ (NR* 1 1 1 ) <FF (I ) * 1 = 1 *80 ) 

1 1 1 FORMAT(80A1 ) 

110 CONTINUE 

Wf?ITE(NW*B8B ) 

noO F ORMAl (//5X* -KCA-Sri Wrix’ IJCK Ml, M*i‘.2 MEfyl Of' CAV*C.33M|.:M PLASX ) 

DO 777 JJKL=l*3 

READ(NR* 101 )A*B*T*EM1 * ETA * DEL.T * UREF * U *RH0 I ,RH02*H 
101 F0RMAT(8F10.5) 

read (NR* 102)1 JK« AKF 
lOa FORMAT! I 1 ♦FlA.a ) 

EMI =EM1 *( .1 )*)t7il>RH01 =RH0l / (12. ) **4$RHO2 = RH02/ (T2. ) **4 

DSTAP=7*/72* J^DELT 

DO 5 L = 1 * 75 

ALl = (L“l )*2.+l • 

no 5 N= 1 *24 

AM1= (N-1 * )-X2, + l • 


AKl = ( AL 1*PI /A )**2+{ ANl *PI/B )*x2 



OMC (Lt N )=CO-»SQRT ( AKl ) 


AK (L *N ) =AK 1 

5 CONTINUE 

DO 6 L= 1 »75» 2 
OMCUT = OMC <1*1 )-5, 

6 AFACT(L)=1.0 
DO 7 ,L=2 *74 *2 

7 AFACT<L)=-1.0 
DO 8 N=1 *24*2 

8 BFACT(N)=1.0 
DO 9 N=2.24*2 

9 BFACT<N)=-1 .0 
WRITE<NW*BOl )U 

801 FORMAT <//20X*Fl5, 6 ) 

DSTl =DSTAR 

Q= 0 , tn-K-RHO 1 *U**2 
C0NST=Q**2*DST1/U*1 • 

CONSTl =0,0064516 
C0N1=DST1/U . 

UC = 0,8*-U 

802 FORMAT ( //20X * -H-LOWER CUTOFF FREO , * ( BULL) * ♦ F 1 5 « 6 * K-UPPER CUTOFFK- ,F 1 5, 
1 6 ) 


W = 0,0 



JK = 0 


D0M=1 .0S0M2 = 2**PI «-DOM 
OMEG=0.O5/CONl 

EEl =EXP (0. 1 ) $EE2 = EXP (0«0235 ) I.EE3 = EXP ( 0 . 40 ) 
DO 2 J=1 «5 
0M3=0M2 
OM2 = OM3 0. 

DO 31=1*90 
JK= JK+1 
AI = I 

OMl=OM3X (1 .^AI/10. ) 

0ML = 0 . I -X0M3 
F=OMl/2./P 1 

OMCl =250. »1 ./CONI *2. *P I 
0MCC = 0MC1 /I 0, 

IF (OMl .GT.OMCC )G0 TO 28 
IF (OMl .GT.OMCl ) GO TO 31 
IF (OMl .LE.OMEG )G0 TO 23 
AT0M = 8« -X-OMl *C0N1 
TW0M=2 .*0M1 *C0N1 
FRACOM = 0.47^^OMl -X-CONl 
I F ( ATOM- 12. )12«13il3 


12 El =1 ./FXP ( ATOM ) 



GO TO 14 


13 E1=0*0 

14 CONTINUE 
IF(TW0M-12. ) 15, l6o 16 

15 E2=l ./EXP(TWOM) 

GO TO 17 

16 E2=0*0 

17 CONTINUE 

IF (FRACOM-1 2 . 118,19,19 

18 E3=l ,/EXP{FRACOM) 

GO TO 20 

19 E3=0.0 

20 CONTINUE 
E=E2*3,7+0.8-XE3-3,4-^E 1 
GO TO 22 

28 E = 3. 1 5/225. 

GO TO 22 
31 E=0.0 

GO TO 22 

23 E= (OMl/OMEG) )f-*4* (3.7/EE1-3.4/EE3 + 0.8/EE2 ) 

22 CONTINUE 

FORE=CONST*E 

803 F0RMAT(//20X»*FORCING FIELD — AMP . =* , F 1 5. 6 , K FREQ * = 1^ i F 1 5. 6 ) 



BK = OM 1 /UCSOK 1 = 1 .+COS (BK*A ) S8K2 = S F N ( BK*A ) 
DU=CMPLX<0.0»0.0) 


ETAl =ETAi>OMl 
OMN= (OMl /CO)**2. 

XX = PI /Aa.XXX = xx**2»ASO = A^<-^<-aSBKSQ = BK-«-K'2 
IFFIJK.EQ.l )PH02=0»0 
DO 10 L=I ,75 
AL= (L-l . )-Ka, + l , 

AL 1 =AL*XX 

APL=A8S ( AL 1 -BK 1 

IF (APL.LEoEPS )G0 TO 51 

ALX1={AL1 «*2-BKS0 )*ASO 

D1 1 =BKl *8, «-AL*ArACT (D/ALX 1 

D1 2=D1 1 ^<bK2/BKl 

GO TO 52 


51 D1 1 =0* 

D12 = -4, /PIX-AFACT(L) 

52 CONTINUE 

DO 1 0 N= 1 ,24 
AN= (N-I * )-«2.+l • 

XI 1 =01 1 /AN«BFACT(N) 
XI 2 = D1 2/AN«-BFACT (N) 
X3=CMPLX (XI 1 ,X12 ) 



AKl =AK (L- N> -OMN 
I F (OMl •L't .OMCUT )GO TO '■Ki 
I F ( ABS ( AK 1 ) .LC *F-PS )GO TO 
1 F ( AK 1 -EPS ) 9 1 1 9^. * 93 

91 AK2=-AK1 

AK2 = S0RT ( Ai<2 ) 

CT=COS ( AK2 ) /S 1 N ( AK2 ) 

XI =- (EMI +RH02*H*CT/AK2 )*OMl •«'*2+AK CL » N )-)^T+RH02-«-G 
XI =X1 +AKF 

X5 = ET A 1 +OM !■«••}( 2-«-RH01*H«-l ,/AK2 
X2=CMPLX (XI «X5 ) 

GO TO 97 

93 AK? = SORT ( AKl ) K-H-X2 * 

CX=EXP( AK2 ) 

CTH= (CXH-1 • )/(CX-l • ) 

XI =- ( (RHOl +RH02-MCTH) *2. >f-H/AK2 + EMl )-«OMl **2 + AK (L ,N T +GKRH02 

XI =X1 +AKF 

X2 = CMPLX (XI »ETA1 ) 

97 CONTINUE 

CMl (L»N)=CABS(X3/X2) 

DU=X3/X2+DU 

92 CONTINUE 
I 0 


CONTINUE 



987 FORMAT ( IbEti. 2) 

W1 C JK ) =CABS ( DU ) **2 

W1 ( JK+1 OOQ )=VM ( JK )-)^FORE 

W1 < JK+1 000 )=CONSTl-»-Wl (JK+1 000 ) 

W=W+W1 (JK+1 000 )*OML 

WR I TE ( NW t 909 ) F , FORE » W I ( JK ) t W 1 (JK+1 000 1 i W 
909 FORMAT <bX.5E20. 7 ) 

3 CONTINUE 
2 CONTINUE 

1 continue 

777 CONTINUE 
STOP 
END 


C GIVE EMI =EM1 *1 0,**7$RH05 =RH0X-1 44 .**2 

DATA TO 8E FED ARE LENGTH OF MODEL ( A 1 , W 1 DTH OF MODEL ( B )♦ TENS I ON OF MEMBRANE(T) 
MASS OF MEMBRANE(EM1 ) «DAMPING(VI SCOUS ) COEFF I C I ENT ( ETA ) » BOUND .L A VER THKNSS(DELT) 
♦REF.VEL (UREF ) 1 VELOC I TY (U) .DENSITY FLD 1 ( RHO 1 ) . DENS .FLD2 ( RH02 > HT OF CAVITY(H) 
FORMAT NUMBER IS 10 1 ♦ 8F 1 0 . A 

26.0 8.0 0.02 1.1 0.0002 0.6 636.0 636.0 

0,00234 1.9379 0.5 


50.0 


18,0 1.0 


0.0002 1.6063 1318.9 1318.9 



U • LHJii 


U f UL'c- 


u » u y 


50.0 iti.O I .0 .1*1 0*0002 

0,00234 0,00234 0,00 

1 0,666666 


1 ,6063 


] 3 1 B . 0 


1318.9 



APPENDIX D 


■ CODE FOR ANALYZING PERIODIC STRUCTURE RESPONSE 



PROGRAM KRAMER ( I NPUT ^OUTPUT * TAPE5= 1 NPUT » T APE6=OUTPUT i TAPE 1 ) 


C f'-x 

C THIS PROGRAM EVALUATES THE MEAN SQUARE RESPONSE OF A KRAMER 

1 PERIODIC SURFACE. THE SUREACE IS A RUUBCR DIAPHRAGM fsUPPORTED 
Z periodically <80 MIL APART) 8Y RUBBER STUBS AND FILLED WITHBWATER. ' 

C-X-* INDUCED aerodynamic EFFECTS USING PISTON THEORY. 

C-it^! EFFECTS OF BACKING USING POTENTIAL FLOW THEORY. 

C** theoretical model USES LINEAR THEORY 

C*X- PROPER values OF DAMPING INCLUDED USING LAB ESTIMATES 

C ANALYSIS SEE SECTION 6. 

C ALL UNITS LB-SLC-FT 

C BULLS SPECTRUM FOR EXCI 1 ING FIELD 

C*-* 

C * -St ***** -it * R #***» -it *•«•**■»****■«■** x- )t * * ■Jt ** )t **** -)t ** -It **■«■*•«■***** X- ** -It -It ***** )t ***** * 

DIMENSION A <21 ) .B (21 ) «C <ai ) . D(21 ) ♦ I J (2! } . FORE ( 1 2000 ) . W ( 1 5 ) » 

1 ASUM (15) 

DIMENSION GE(2I ) ,F(21 ) 
dimension HEF(21 ) 

NR=5$NW=6 
C DATA ■ 

H=5O,0A0$RHO=62.4SG=386*giU=720.3>ETA = 0. 00 1 3.E = 2000.$EM1 =8 .8$ANU=0. 45 



DELTA=1 .08SY = 0,3>TP I =6 • 283 1 B5260-‘liCO = 57000 , 

ETA=0. 0000001 
EPS=0. 0000000 1 

C INITIALISATION OF DATA 

EMI =EM1*( ,1 )**63>RHO = RHO/32*2/12./12«/12./12.3;DST1=7*/72.*OELTA 
AKF = E»A i ==H/AKF$Pl = A1 
AKF = 20.03.Al=H/AKF3iPl =A1 
EMI =EM1 *A I 

DD = E*H«'-k-3/3./ ( I .-ANU**2) 

D1 =DD*A1SETA1 =ETA*A1 1UC = 0 ,8*U3>Q = 0«5*RH0*U**23>C0N1 =DST1/U$01 =Q*1(-2-)(C 
1 ONI 

G1 =G*RH0*A1 

OMEG = 0 . 05/C0N1 $D0M=0 . 1 -X-TP I aiOM2 = DOM 

D0M=1 O.-X-TPI 

0M2=D0M 

OMCl =25.*TP1 *U/DELTA 
0MCC=1 0.*OMC1 
C initialise ARRAYS 

DO 101 J=l,15 

ASUM ( J ) =0iW ( J )=0, 

1 01 CONTINUE 

C MAJOR LOOP for increment OF OM 


JJ = Oi&EE 1 =EXP <0. 10 )iEE2 = EXP ( 0.0238 >iEE3 = EXP (0.4 0 ) 



C DO 2 JK=1 <6 

OM3=OM2aOM2=OM3*l 0, 

C DO 1 JL=1 i90 

DO 1 JL=l ,20000 
AJsJL 

OMl =0M3^e ( 1 .+AJ/10. ) 
FF=OMl /TP I 
JJ = JJ-H 

uc=o •a«-u 

IF (OMl «UE.OME:.G)GO TO 23 
IF (OMt .GE.OMCl )GO TO 5 I 1 
IF (OMl .GT.OMCOGO TO 513 
AT0M=6. 0*0M] *C0N1 
TW0M=2**0Ml *C0N1 
FRAC0M=0.47*0M1 -K-CONl 
IF (ATOM-60. ) 1 2, 13 , 1 3 

12 El =1 ♦/EXP (ATOM) 

GO TO 14 

13 E1=0. 

14 CONTINUE 

1 F (TWOM-60^ ) 1 5, 16 1 1 6 

15 E2=l ./EXP(TWOM) 


GO TO 17 



16 E2=0.0 

17 CONTINUE 

IF (FRACOM-60. ) 18, 19i 19 
la E3=l »/EXP(FRACOM) 

GO TO 20 

19 E3=0. 

20 CONTINUE 
E=F2*3»7+0.8*E3-3«4*E1 
GO TO 22 

511 E=3,1B/225. 

UC=0»39itU 
GO TO 22 
513 STOP 

23 E=3,7/EE1 +0.a/EE2-3*4/EE3 

E=E* (OMl/OMEG )**4 
22 CONTINUE 

FORE (3J )=Q1*E*0. 0000 1 

AKl = (OMl /CO )**21>AMU'=0M1 /UC^|-I*ETA I =ETA 1 *0M 1 
M=-l 1$SUM=0.£SUM1=0. 

DO 30 Ml=l ,21 
M=M1 -1 1 
AM 1 =M 


AML= ( AMU+ApM 1 -X-TP I ) /H 



AK = AML**2-AK1- 


A{M1 ) = Dl*AML**A + Gl“OMH(-)^2)tEMJ3:B(Ml )=ETAF 
IF<ABS(AK ).LE*EPS)G0 TO 40 
lF(AK-EPS)39.40i41 
40 IJ<M1 ) = 1 
GO TO 47 
39 CONTINUE 
I J (Ml ) = 0 
AK=-AK 
AK=S0RT < AK) 

,SK = AK^H 

SI =SIN (SK )SC1 =COS(SK) 

IFFSl.EQ.O) GO TO 901 

A(Mi )=A(Mi )+Rh 0*A1*C1/SI/AK*0M1**2$B (Ml )=B(Ml ) +RhO*A 1 /AK*OMl 
GO TO 902 

901 I J (Ml ) = 1 

902 continue 
GO TO 47 

4 1 I J (M 1 ) =0 

AKsSQRT ( AK ) 

% 

SK=AK*2 .XH 
SK1=EXP (SK) 


CTH= (SKI +1 . )/ (SKI -1 . ) 



52 CONTINUE 

A<M( )=A (Ml 1 -OMl**2-«'RHO/AK*Al * ( 1 .+CTH) 

47 CONTINUE 

IF ( F J(M1 ) )48i48«49 

48 CONTINUE 

A2 = A(M1 )*X-2+B(Ml ) **2 
SUM=SUM+A (Ml )/A2 
SUMl =SUM1+B (Ml )/A2 

49 CONTINUE 
30 CONTINUE 

A9 = A ( 1 1 ) )t)t2H-0 ( 1 1 

GAMl = A I » ( A { 1 1 )* (SUM+ 1 • ) +f3( 1 1 )-XSuM 1 ) 

GAM! =GAMI /A9/ ( { 1 . +SUM ) * X2 + SUM 1 ) 

GAM2 = A1/A9* ( A ( 11 )*SUM1 -B (11 ) ^( ( SUM+l . ) } / (SUMl **2+ C 1 .+SUM>* X2 ) 
DO 60 Ml=l .21 
IF(IJ(M1 )161 .62.61 

61 CONTINUE 
D ( Ml 1=0.0 
C(M! 1=0.0 
GO TO 67 

62 CONTINUE 

C ( Ml ) = ( ft (Ml )XGAMi +B ( Ml )*gaM 2 ) / i a (M] ) ■*■*£; + tM Ml ) « > 

D (Mt ) = ( B (Ml ) X-GAMl -A ( Ml ) -«-gAM 2 >y ( A (Ml ) X*2+B ( Ml ) ) 



92 CONTINUE 

WR ITE (NW»200 )FF> ASUM (A ) «GE U 1 ) ♦ GE ( 1 2 ) ♦ GE ( 1 0 ) »HEF (1 1 ) « HEF M2) «HEF( 1 
10) ,W(4) ♦FORE(JJ) 

603 F0RMAT</50X4^< IM* ) 

209 FORMAT ( 1 XM 0E12 .5 ) 

208 FORMAT (//5X. 10E12 .6 ) 

1 CONTINUE 

2 CONTINUE 

201 FORMAT! 1X»7E16.T) 

STOP 

END 



APPENDIX E 


CODE FOR ANALYZING MEMBRANE UNDER NARROW AIR GAP 



JOB ♦ 1 4 £500 4 77000 4 30000 


A4677 


R4 6 SB 


1 0071 8 


B I N34 


RUN (S ) 

REQUEST 4 TAPE 1 4 HY , SA VTP 4 R I L 4 RS8 4 TW2 0 . 

REWINn(TAPEl ) 

PROGRAM MEMGAP ( INPUTt OUTPUT 4 TAPES= I NPUT 4 TAPE6 = 0UTPUT 4 TAPE 1 ) 
DIMENSION W( 51 4 51 ) 4 WO <51 451 ) 4 WN < f5 1 4 5 1 ) 4 P ( 5 I 4 5 1 ) «AU( 10! ) 

C T= TENS ION (LBF/IN) 

C RHO=MEMBRANE DENSITY (LBM/CU.FT) 

C E=SU0STRATE MODULUS OF ELASTICITY (PSIl 

C D=S<mSTRATE DEPTH (IN) 

C GAP=GAP BETWEEN MEMBRANE AND SUBSTRATE (IN) 

C H=MEMBRANE THICKNESS (IN) 

C B=DAMPING COEFFICIENT ( L 8 F-SEC/CU . FT ) 

C DX=ELEMENT SIZE (IN) 

C NI=NUMBER OF NODES IN X-D I RECT I ON 

C NJ = NUMBER OF NODES I N, Y-D I RECT I ON 

C NT=NUMBER OF TIME STEPS — nT = nx/C WHERE C IS WAVF-' SPE-TD 

C GM IS TIME STEP CONTROL DT = DX/ ( C-i< CM ) 

T= 1 

RH0=87 
F= 1 . 





NW=6 


H^O.OOl 

B=o*oooobi 

N I =B] ;r.NJ = '51 ‘f-DX = 0.2 
U = SO't.f)FLTA= 1 4l«Pl=3* 

NT=2n000 

PMAX = 0*00234*U*^(-2*0.5*0,01 
P) =PMAX-KU**2/DX*!(2 
PI =p 1 * I 44 . 

F= 1 .?*UK-12./PI/2./nELTA 
C convert INPUT TO CONSISTENT UNITS 

T=12.*T 

RH0=RH0/32. 1 74 
E=144.«-E 
D=D/ 1 2 
GAP=GAP/12. 

H=H/ 1 2 . 

C WAVE SPEED IS GIVEN BY SO. RT. OF T/(RHO«H) 

C= (T/ (RHO-KH } )**0,5 
OX = DX/ I 2 
DT = SHO » 1 *-)f6 


gm=ox/c/dt 



I \M T I L ! 7"“ D I SPl, ACF N'€NTF 


DO 7 ! - 1 1 (■! I 

DO ) J= 1 * NO 
WO ( I , J )=0 . 
l>/ ( r 1 J ) = 0 . 

M . J ) =0 , 

P ( 1 . J ) -• o . 

COMT 1 NUF 
CONT i NUf 
B = R r 

CALCnuATF NFVJ D [ SPL ACF.’MPNTS 

PRFF‘ loading can bi. fntff^ed here 

N I y~N ! - ] 

NJY=NJ- 1 

A 1 =B*GM/ ( 2 , -Jf-T-Jf-DX ) + { GM/DX ) **2 
A2= 1 ./DX^( 1 ./OX 

A3 = fi>iGM/(2. imiDX )- (GM/nX) If- !( 2 
AA=F/(r)-«T 1 
AE= 1 ./T 

TF= I • /F.'i.OM = ? . -)f-PI >fF3.T 1 =OX/U-'tTWOT I =2 , H-T 1 
TPLTl =TE+T] 


TPLT2=7 F + TWOTl 



ENI=NI-1 


?.?1 

£!0l 


1 54 
1 55 

221 

151 
1 56 


rMAx=iiMi/‘- "|■'X+l v/r-ifix/u 

REWIND 1 

DO 1 30 N=1 iNT 

I T=2+N 

T I mf-NX-D I 

T=T IME 

no PF7 !X=-P,N'I 

DO ??7 JY-?,NJ 

P ( I Xi JY )--n.o 

I F (TI MF-T-.1AX )?.0 1 ,210,210 

DO 209 IX--2,N!X 

EXX= 1 X-2 

T=T IMF--r XX!<|1X/U 

I F ( T ) 1 52 , t ■■•'2 1 I 54 

1F<T“T1 ) ! 55i 155* 151 

CONTINUE 

DO 221 JY=2,NJY 

P( IX, JY )=Pl->f (T/0M~SlN(0M)fT)/0M««2) 
GO TO 152 

IF'< T-TWOT 1 ) 1 56, 1 56, 1 57 
CONTINUE 


DO 222 JY=2,MJY 



?22 P(rX*JT')=Pl«< TWOT 1 / 0 <v-T/ 0 MH-? 2 . XS 1 N (OM« <T-T I ) ) /OM-X--X-2-S I \i (OMK T ) /om K-2 
1 ) 

GO Tf) 1 

157 1 F ( T-TE ) 1 , 1 58 1 1 50 

158 CONTINUE 

DO 223 JY=2,NJY 

223 P ( r X » JY ) = Pl «2 * /OM-**2-«'S 1 N { OM* (T“T 1 )) * ( 1 . -COS < OI^Ht T 1 ) ) 

GO TO 1 52 

150 IF ( T-TPLT 1 ) 1 80 , 1 60 . 1 81 

160 CONTlNUr- 

DO 224. JY=2.NJY 

224 P ( IX , JY )=Pl iM < T-Tl ) /OM+2 * /0M**2*S IN ( 0M«- ( T-T 1 •) )~S IN ( OM X' ( T- TWOT 1 ) )/0 
1 M* X? ) 

GO TO t 52 

161 IF (T-TPLT2 ) 162, 162 ♦ 1 52 

162 CONTINUE 

DO 225 JY=2,NJY 

225 P( IXi JY)=P1 X( (T-TE-Tv;0T1 ) /OM-S I N ( OM* ( T-TWOT 1 ) )/0M**2 ) 

152 CONTINUE . 

209 CONTINUE 

210 CONTINUI- 

DO 20 I =2, NIX 


DO lO J=2»NJY 



SUM = '*' ( 1 + 1 4 J )+W ( I“1 ♦ J )+W ( I ^ J-l- 1 )+\> { I 4 J-1 ) 


SUM='^UM+ (2 • *GM*#2-4 . )*W < ! 4 J ) 

SUMj =W0 { I ; J > 

SUM2 = A4*(U'( 1 4J)+GAP) 

WN M 4 J ) =SUM*A2/A1 +SUMI-X A^/Al ~P ( I 4 J )*A5/A 1 
I F ( W < I 4 J )+GAP>S 4 5 46 
i"N( I 4 J ) =2.*W0 ( 1 4 J )-W ( I 4 J ) 

CONT I NUE 
CONT I NUF 
CONT I NUF 
A I =N1 - 1 
BI=] 00. 

DEL=BI/A1 
1DFL=DEL +0*2 
DO 19 1=1,101 

El = I-l 

J= ( I-l l/IDEL+l 
EJ=J- 1 
J1 = J+ I 

X I =F-I /DEL-EJ 

AU( I )=WM( J, J 1 )-«•( } ,_xi )+WN( J1 4 n ) XXI 
CO^lT I NUF 


WR I TF ( 1 ) AU 



c 


SHIFT LOCATIONS OF W ARRAYS 


DO 91 I=?iNIX 
fv. . )i J = ;-',NJY 
'.HO ( I ♦ J ) = l»' ( I ♦ J ) 
Ui ( I . J ) " WN ( I 1 J 1 
91 CONTINUE 
130 CONTINUE 
STOP 


END 
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